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Abstract 

We apply the generalization of the Parikh-Wilczek method to the tunneling of massive 
particles from noncommutative inspired Schwarzschild black holes. By deriving the equation 
of radial motion of the tunneling particle directly, we calculate the emission rate which is 
shown to be dependent on the noncommutative parameter besides the energy and mass of 
the tunneling particle. After equating the emission rate to the Boltzmann factor, we obtain 
the modified Hawking temperature which relates to the noncommutativity and recovers the 
standard Hawking temperature in the commutative limit. We also discuss the entropy of 
the noncommutative inspired Schwarzschild black hole and its difference after and before a 
massive particle's emission. 
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1 Introduction 



Black holes are an incredible outcome of the Einstein's theory of gravity. They were thought 
that no matter inside could escape and thus they were invisible from outside. In 1970's, 
Hawking startled the physical community by proving [1] that the black holes are not really 
black. They can radiate like a black body with the "Hawking temperature", T H = ^, where 
k is the surface gravity of black holes. Although Hawking first suggested the heuristic picture 
of the radiation as tunneling, his calculation was completely based on quantum field theory 
in the curved spacetime, which is independent of a tunneling process. 

An intuitional method was proposed at the end of 1990's by Parikh and Wilczek [2] in 
which the radiation is dealt with as tunneling of particles. It is suggested in the method that 
the barrier is created by the tunneling particle itself. With the WKB approximation, the 
emission rate is calculated semiclassically. In terms of the equality of the emission rate and 
the Boltzmann factor, the Hawking temperature is then derived. The results show that the 
emission spectrum is not precisely thermal but its leading order coincides with the standard 
Hawking radiation, i.e. the black body spectrum. However, it is only investigated in ref. [2] 
that massless and uncharged particles tunnel from spherically symmetric black holes, such as 
the Schwarzschild black hole and Reissner-Nordstrom black hole. When the tunneling of a 
massive and charged particle is considered, the problem that the trajectory of the tunneling 
particle is not null geodesic has to be confronted and thus the equation of motion of the 
massive charged particle has to be found out. In ref. [3], the concepts of phase velocity and 
group velocity are introduced in order to derive the equation of radial motion for a massive and 
charged tunneling particle, and the equation of radial motion is obtained through equating 
the radial velocity of the tunneling particle to the phase velocity. In addition, a quite different 
method based on the Hamilton- Jacobi equation is provided in refs. [4, 5, 6, 7]. This method 
is not involved in the trajectory problem, and can thus be applied to tunneling particles with 
mass and/or charge naturally. Recently, the Parikh- Wilczek method has been generalized [8] 
to the tunneling particle with mass and charge in terms of strictly deriving its equation of 
radial motion. The corresponding results show that the emission rate depends on the mass 
and charge of the tunneling particle, and that the radiation temperature relates to the charge 
but equals the standard Hawking temperature in the limit of zero charge. 

On the other hand, the idea of noncommutative (NC) spacetimes has been revived since 
the publication of the Seiberg and Witten's work [9]. This idea is originally expected to 
remove the divergence in quantum field theory [10] and later expected to be an indispensable 
ingredient for quantization of gravity [11]. Till now a large amount of papers have been 
devoted to the construction of quantum field theories on noncommutative spacetimes, for a 
review, see ref. [12]. Gravity theories on noncommutative spacetimes have extensively been 
studied, for instance, see the review article [13], and black hole solutions have also been 
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obtained [14]. Most of the works related to NC spacetimes are based on the "star-product 
approach" in which the noncommutativity of spacetimes is encoded through an ordinary 
product in the noncommutative C*-algebra of Weyl operators to a noncommutative star- 
product of functions in commutative C*-algebra. Within this framework one calculates only 
order by order in noncommutative parameters and therefore loses the nonlocality of non- 
commutative theories. Nevertheless, the so-called "coordinate coherent states approach" has 
been proposed [15] from a different point of view for the study of noncommutative quantum 
mechanics and quantum field theory. The models established in quantum field theory with 
such an approach are consistent with the Lorentz invariance, unitarity and UV-finiteness. 
The main idea of the approach is that the point matter distribution should be depicted by 
a Gaussian function rather than a Dirac delta function. In terms of this smeared matter 
distribution, the NC inspired Schwarzschild black hole solution is provided [16] under the 
assumption that the Einstein's equation is not influenced by the noncommutativity. The 
remarkable property of the solution is that there exits a minimal mass M under which no 
horizon can form and therefore there will be a remnant after the Hawking radiation, which 
may probably solve the "information loss paradox". Based on ref. [16], a lot of efforts have 
been devoted to extend the solution to higher dimensions and to charged and rotating black 
holes, such as for the NC inspired Schwarzschild black hole solution in higher dimensions [17], 
for the NC inspired Reissner-Nordstrom black hole solution [18, 19], and for the NC inspired 
Kerr(-Newmann) black hole solution [20, 21]. Furthermore, the thermodynamics of the var- 
ious NC inspired black holes mentioned above has been analyzed in accordance with the 
Parikh-Wilczek method or Hamilton- Jacobi method, in which only the tunneling of massless 
particles is considered [22]. In addition, the massive particle's tunneling from NC inspired 
black holes is discussed in ref. [23] in terms of the method given by ref. [3]. 

Our motivation is obvious from the above summary on both the aspect of the Hawking 
radiation as a tunneling process and the aspect of the tunneling from NC inspired black holes, 
that is, we follow our generalization of the Parikh-Wilczek method [8] to study the tunneling 
process of massive particles from the NC inspired Schwarzschild black hole. In the next 
section, we give a brief introduction to the NC inspired Schwarzschild black hole. In section 
3, starting from the Lagrangian of the tunneling particle, we derive its equation of radial 
motion directly and then calculate its emission rate. After equating the emission rate to the 
Boltzmann factor, we therefore obtain the modified Hawking temperature which relates to 
the noncommutativity of spacetimes. In the following section we discuss the entropy of the 
NC inspired Schwarzschild black hole and its difference after and before a massive particle's 
emission. Finally, we make a conclusion in section 5. 
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2 Noncommutative inspired Schwarzschild black hole 



Consider a noncommutative spacetime where 9 denotes the noncommutative parameter. In 
the coordinate coherent states approach [15], a particle (black hole) with mass M is smeared 
over a region with width \/9 and is described by a Gaussian function rather than a Dirac 
delta function as follows, 

M _r£ 

Pe(r) = w e 49 • (!) 

(4tt#)2 

If the Einstein's equation is not influenced, the Schwarzschild-like solution of the Einstein's 
equation with the matter source mentioned above takes the form [16], 

ds * = - (i - 2 -Ml\ & + (l - ir > + rW, (2) 

where the parameter Mg with the dimension of mass satisfies the following formula, 

Pe (r')Anr' 2 dr' = — 7 (^J " ( 3 ) 

The gamma function in the above equation is defined to be 7(0, 6) = y t a e~ l . 

Here we make a comparison between the coordinate coherent state approach and the 
star-product approach. For the former which is argued [24, 22] to be equivalent to the so- 
called Voros- product (a specific star-product), the effect of noncommutativity only on the 
matter is considered while the gravity is left untack. Such a treatment results in a smeared 
"point mass" described by a Gaussian function (see eq. (1)), but it maintains the usual 
Einstein equation unchanged. For the latter, the noncommutative effect on both the matter 
and gravity is considered, and therefore it is in general hard to find exact solutions to the 
noncommutative Einstein equation. Only in some very special cases can the nonperturbative 
solutions be found, such as that discussed in ref. [25] (mentioned by ref. [26]). Based on 
the twisted noncommutative gravity theory [27], the nonperturbative solutions are derived 
under some specific compatibility conditions of metrics and twists. There is no proof that 
the coordinate coherent state approach yields a solution of the noncommutative Einstein 
equations, such as that proposed in ref. [27] which works under the star-product approach. 
Comparatively, it is not so hard to find exact solutions and to study the thermodynamics of 
various noncommutative inspired black holes under the framework of the coordinate coherent 
state approach. In some sense, the coordinate coherent state approach might be regarded as 
an effective way at least from a phenomenological viewpoint. As a whole, the both approaches 
corresponding to different star-products have their own merits in light of the present studies 
and more researches on various black holes are needed if one tries to distinguish a better one 
from the two. 
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The solution (eq. (2)) describes a spacetime interpolating a de Sitter spacetime near the 
origin and a Schwarzschild spacetime far away from the origin. Taking the limit r — > oo in 
the gamma function, one can see that the parameter Mq is just the total mass M of the 
source. In this limit, the solution (eq. (2)) coincides with the classical vacuum Schwarzschild 
solution. Alternatively, we can get this result by taking the limit 9 — > 0, which presents the 
degeneration from the noncommutative to the commutative case. Furthermore, the square 
of line element has no singularity at the origin, which can be shown by the Ricci scalar near 
the origin, 

AM 

R(0) = "=s ■ (4) 

If the equation goo(r H ) = has solutions, there exist horizons. That is, the horizon r H is 
determined by the following relation, 

, N AM /3 r%\ 

However, it cannot be solved analytically. The numerical analysis shows [16] that only when 
the total mass M is larger than the minimal value, M ~ 1.9y/9, can the horizons be formed § . 
This means that there will be a remnant after the end of Hawking radiation. 

Due to the fact that 7 function tends to unit quickly, the ordinary Schwarzschild horizon 
r' H = 2M is a good approximation. To have a more precise approximation to reflect the 
noncommutative effect, one substitutes r' H = 2M into eq. (5) and thus has the relation, 

/ / M\ 2M m 2 \ 
r H = 2M e (r' H ) = 2M (ert ^ j - , (6) 

where the Gaussian error function is defined as erf(x) = ^= J* e - * 2 ' dt. When we repeat this 
procedure, th will get closer and closer to the true value. In this paper we only consider the 
first iterative approximation given by eq. (6). 

3 Quantum tunneling of massive particles 

Now we turn to analyze the tunneling process of massive particles by following the idea of the 
Parikh- Wilczek method [2] together with its generalization [8] . The first thing we have to do 
is to derive the equation of radial motion of the tunneling particle. As we know, when a free 
massive particle with mass m moves along a time-like geodesies, its trajectory is determined 
by the Lagrangian, 

dx^ dx u 

2C = mg, u — — , (7) 



§ The authors of ref. [28] use modified dispersion relation to analyze the thermodynamics of the noncom- 
mutative inspired Schwarzschild black hole and give some interesting results which show that the remnant 
mass may be larger than M . 
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where the parameter r is the proper time. The choice of this parameter provides an extra 
equation, 

g^x" = -1, (8) 

where a dot stands for the derivative with respect to r. Substituting the metric (eq. (2)) into 
the Lagrangian (eq. (7)), we have 

From the Lagrangian we know that the canonical coordinate t is a cyclic coordinate. As a 
result, the corresponding canonical momentum, i.e. the particle's energy is conserved, 

dC / 2M e \ . 
— p t = ~ ~7^~ — mil I t = u = const., (10) 

where the minus sign before p t presents due to the positivity of the energy ui of the tunneling 
particle. Substituting eq. (10) into eq. (8), we obtain 



r = ±m- 1 ^cu 2 -m 2 (l-^^J , (11) 

where the upper (lower) sign in eq. (11) corresponds to an outgoing (ingoing) geodesies. 
Using eqs. (10) and (11) we finally derive the equation of time-like geodesies outgoing along 
the radial direction, 



* '' '^^L.rf/,-*, (12) 



dt t uj \ r 

In order to investigate the Hawking radiation as a tunneling process, we have to transform 
the coordinates into the ones which have a good behavior at the horizon. As done in ref . [2] , we 
choose the Painleve coordinates [29] which can be acquired by the following transformation, 



dt p = dt+ v/1 900 dr, (13) 

#00 

where t p denotes the Painleve time and the other components of coordinates keep unchanged. 
Under this coordinate transformation, the square of line element (eq. (2)) becomes 



ds 2 = - (l - dt 2 p + dr 2 + 2^^dt p dr + r 2 dQ 2 . (14) 

Correspondingly, the equation of time-like geodesies (eq. (12)) in the Painleve coordinates 
takes the form, 

-l 



dt. 




^ - m? (1 - !?■) V r 



In accordance with the WKB method, we have to get the imaginary part of the action of 
the tunneling particle for the calculation of the tunneling rate. This imaginary part emerges 
from the integral of the radial coordinate if the corresponding integrand has a singularity, 
and it is defined [2] for the tunneling particle which crosses the horizon outwards from r- m to 
r out as follows, 



ImS 



Im 
— Im 



rr ut /-rout rp r pr ou t r 

/ p r dr = Im / / dp r dr = Im / / 

J r-,„ J r\„ JO J V;„ J rr 



dH 

dr 

m dt v 



dr 



rout /•« d ~ 



dr 
m dt p 

dH 



(16) 



where the Hamilton equation JF- = |~- has been utilized, and the minus sign appears due to 
the relation [2]: H = M — Co. As explained in ref. [8], the lower bound of the energy is m 
instead of zero because the tunneling particle is massive, which leads to a modification to the 
emission rate. When we consider the self-gravitation of the tunneling particle, we have to 
make the replacement M$(M) — > Mg(M — u) in eq. (6) and eq. (16). Moreover, with eq. (5) 
we give the initial and final positions of the tunneling particle as follows, 



4M 



7T 



3 r 



2 
in 

2' 46 



Tout = 



4(M-u;)_ (3 

2' AO 



71 



-7U 



(17) 
(18) 



At present we exchange the order of integration in eq. (16) and pick out the r integral, 

-i 



Im 



= Im 



dr 
dt n 



dr 



2M e (M - u) 



+ 



uS 1 — m 2 



1 - 



2M b {M-Cj) 



2Mb (M - u) 



dr. 



J 



There is a pole at r = 2M(M — Co) in the integrand. Integrating r by deforming the contour 
and making use of the residue method and eq. (6), we obtain 

-l 



Im 



dr 
dt n 



dr = -AttM 9 (M-Cj). 



(19) 



Next we complete the energy integral in eq. (16) and at last work out the imaginary part of 
the action, 



ImS = 47r M e (M-u)duj 



(M-uf--e 



-2tt 
+2vr 



{M-mf--9 



erf 
erf 



M — U \ / (M-u) 

■ - 6v^ (M - u) e'—^ 



V9 



(20) 
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The tunneling probability takes the form, 

r~e- 2Im5 ~e ASBH , (21) 

where ASbh is the difference of black hole entropy after and before the Hawking radiation. 
From eqs. (20) and (21) one can see that the emission rate depends not only on the energy of 
the tunneling particle but also on its mass m and the noncommutative parameter 9, and the 
dependence of 9 will lead to a modified Hawking temperature. Expanding ImS with respect 
to ui and m to their first orders, we get 



ImS = 4ttM 



JM\ 2M _Mi 

erf — — e « 

\V9j V^9 



(u-m) + 0(u,m). (22) 



Equating the tunneling probability T ~ e 2ImS to the Boltzmann factor e t h , we can deduce 
the modified Hawking temperature from eqs. (21) and (22), 

1 



8nM 



erf 



M \ _ 2M_ 



(23) 



which is valid only when M is larger than M , the mass of the extreme black hole. When 
substituting eq. (6) into the above equation, we obtain a simple relation between the modified 
temperature and the event horizon: T H = which takes the same form as that of 

the ordinary Schwarzschild black hole. However, we note that r# is just the first iterative 
approximation for the NC inspired Schwarzschild black hole while the same symbol stands for 
the exact value of event horizon for the ordinary Schwarzschild black hole. In the commutative 
limit 9 — > 0, eq. (23) recovers the standard Hawking temperature, T H = g^g, for the ordinary 
Schwarzschild black hole. 



4 Entropy of noncommutative inspired Schwarzschild 
black hole 



In this section, we discuss the entropy of the NC inspired Schwarzschild black hole and its 
difference after and before a massive particle's emission. In the case of massless tunneling [2], 
ImS = — | ASbh- But in our case, this equality no longer holds. We calculate the entropy of 
the NC inspired Schwarzschild black hole by using the first law of thermodynamics, dM = 
TdS, 

" M dM' 



= / 

J A 



M 

An 

-47T 



Ti 



H 



M 2 --9 


erf ^ 


2 _ 





Ml 



2 



M 

-3 



12Vtt0 Me 



12Vn9M Q e- — , 



(24) 
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where the temperature given by eq. (23) has been used. Note that the lower integration 
limit cannot be set to zero but to the mass of the extreme black hole Mo because the NC 
inspired Schwarzschild black hole can exist only for M > M . This is a basic feature which 
we should consider in the calculation of the entropy, however, it was unnoticed by some 
literature [22, 23]. The entropy can reduce to the standard Bekenstein-Hawking entropy, 
>Sbh = AnM 2 , for the ordinary Schwarzschild black hole in the commutative limit 8—^0, 
in which case no extreme black holes exist. With eq. (24) we obtain the difference of the 
entropy after and before the emission, 



ASbh = S bh (M-uj) 



(M-u) 2 --8 



Sbh(M) 
erf 



M 



V8 



12Vn8 (M - uj)e ~ 



An 


M 2 --8 


erf ^ 




2 _ 





erf i-j=\ - nV^Me" — , 



(25) 



which is independent of the mass of the extreme black hole. Comparing it with eq. (20), 
we see ImS' ^ — ^ASbh- When we take the limit m — > 0, which corresponds to the case of 
massless particle's tunneling, the tunneling probability equals the exponent of the difference 
of entropy, e~ 21mS = e ASBH . As a consequence, our result is consistent with that of ref. [2]. 



5 Conclusion 

In this paper, we analyze the tunneling process of massive particles from the NC inspired 
Schwarzschild black hole following the way given in our previous work [8] which is based on the 
Parikh-Wilczek method [2]. We provide some novel and previously unnoticed properties for 
the noncommutative inspired Schwarzschild black hole, such as: (i) We generalize the Parikh- 
Wilczek method to massive tunneling particles in methodology; (ii) The tunneling rate of 
massive particles is suppressed due to the mass, and the effective radiation temperature is 
modified by the noncommutativity but recovers the standard Hawking temperature in the 
commutative limit, and eqs. (20) and (21) imply that the radiation is not precisely thermal, 
as pointed out in ref. [2]; (iii) The well-known conclusion that the tunneling rate equals 
the exponent of the difference of the entropy does not hold in the massive case but it can 
be recovered in the zero mass limit, which coincides with that of ref. [2]. This inequality 
may imply that there exist correlations between tunneling particles. To see this, as done in 
refs. [22, 23], we can compare the probability of tunneling of two particles of energies (masses) 
Co>i (mi) and oo 2 (717-2) with the probability of tunneling of one particle with the energy (mass) 
summation uii + ui 2 (mi + m 2 ). We obtain an inequality, 

In ^ui,mi hi Tui2,m2 7^ ^^-^ u)i+uJ2,mi+m2 i (26) 
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for the massive case, but an equality, lnr wi + InT^ = lnr wi+W2 , for the massless case 
mi = m2 = 0. Therefore, the radiation may also contain information beyond the remnant 
for the massive case, which deserves further investigations. 
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